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Abstract 

We perform total energy calculations based on the tight-binding Hamiltonian 
scheme (i) to study the structural properties and energetics of the extended 
{311} defects depending upon their dimensions and interstitial concentrations 
and (ii) to find possible mechanisms of interstitial capture by and release from 
the {311} defects. The generalized orbital-based linear-scaling method imple- 
mented on Cray-T3D is used for supercell calculations of large scale systems 
containing more than 1000 Si atoms. We investigate the {311} defects sys- 
tematically from few-interstitial clusters to planar defects. For a given defect 
configuration, constant temperature MD simulations are performed at 300 - 
600 K for about 1 psec to avoid trapping in the local minima of the atomic 
structures with small energy barriers. We find that interstitial chain struc- 
tures along the (011) direction are stable interstitial defects with respect to 
isolated interstitials. The interstitial chains provide basic building blocks of 
the extended {311} defects, i.e., the extended {311}defects are formed by 
condensation of the interstitial chains side by side in the (233) direction. We 
find successive rotations of pairs of atoms in the {011} plane are mechanisms 
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with a relatively small energy barrier for propagation of interstitial chains. 
These mechanisms together with the interstitial chain structure can explain 
the growth of the {311} defects and related structures such as V-shape bend 
structures and atomic steps observed in transmission electron microscopy im- 
ages. 

61.72.Ji, 61.72.Nn, 71.15.Nc, 71.15.Fv 
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I. INTRODUCTION 



The decreasing size of semiconductor devices requires precise control of device structures, 
particularly dopant distributions. Ion implantation introduces energetic charged atomic 
particles in a substrate for the purpose of changing electrical, metallurgical and chemical 
properties of the substrate. The wide use of ion implantation is due to its precise control 
over total dopant doses, depth profiles and the area uniformity. However, ion implantation 
induces transient enhanced diffusion (TED) of dopants: the diffusivity of the dopants is 
abnormally enhanced for a transient period of time after ion-implantation.0 Many studies 
have suggested that lattice damages, introduced during ion implantation, are responsible for 
fast diffusion of dopantsBi 

In particular, the diffusivity of boron during annealing in ion-implanted samples is en- 
hanced by many orders of magnitude greater than B diffusivity in a thermal equilibrium.il 
This transient enhanced diffusion of B places limitations of the use of B + -implantation in 
fabricating submicron devices. Boron TED occurs due to excessive Si interstitials, which 
are created by B + -implantation and contribute to the fast diffusion of B.Bi The boron dif- 
fusion via pairing of B and Si interstitials has been supported by experiments and theories.0 
Experiments observe enhancement of the B diffusivity when Si interstitials are selectively 
injected by surface oxidation! and by Si + -implantation.i First principle calculations reported 
that the activation energy of B diffusion associated with Si interstitials is the lowest among 
possible diffusion mechanisms.!'! Therefore, the concentration of Si interstitials after the 
ion implantation and during the thermal annealing process is an important parameter to 
determine the final depth profile of B in the ion-implanted samples. 

The concentration of Si interstitials depends on (1) the ion implantation conditions - 
the incident energy and the dose - and (2) the temperature at which thermal annealing 
is performed. Implantation at a higher incident energy or with a higher dose of the ions 
generates more Si interstitials. As the name transient enhanced diffusion indicates, the 
enhancement factor of the B diffusivity decays as the thermal annealing proceeds, and 
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eventually the diffusivity of the B converges to that in a thermal equilibrium. The temporal 
extent of the enhanced diffusion of B decrease at higher annealing temperatures, since the 
excessive Si interstitials migrate to the regular lattice sites in a shorter period of time.i It is 
observed that the temporal extent (duration of the boron TED) is exponentially activated. 
The activation energy, depending on the ion implantation conditions, shows a large variance: 
1-5 eV.ii 

This activation energy summarizes complex processes for B diffusion in ion-implanted 
samples - generation, diffusion and annihilation of the Si interstitials as well as the inter- 
actions of B with the B and Si interstitials. Recently, Eaglesham et a/.0 and Stolk et 
suggested that emission of Si interstitials from particular extended defects, namely {311} 
defects, causes the boron TED. They observed that the dissolution of the {311} defects oc- 
curs at the same time and the temperature conditions as boron TED. The extended {311} 
defects are detected under B + -implantation at an incident energy of few tens of and 
the corresponding activation energy of the boron TED was obtained at 3-5 eV.@ The high 
value of the activation energy can be related to the stability of the {311} defects with respect 
to isolated interstitials. However, boron TED has been also observed without detecting any 
macroscopic defect, when the energy of the B + -implantation decreases below 10 keV.I A 
lower activation energy (about 1 eV) was estimated for this low-energy ion implantation. 
This low-energy B + -implantation experiment suggests that the formation of stable {311} 
defects depends on the interstitial concentrations. Furthermore, the observation of boron 
TED in the absence of the observable {311} defects suggests that microscopic interstitial 
clusters may exist and contribute to the B diffusion by releasing interstitials at a lower en- 
ergy cost. It is possible that the microscopic interstitial clusters are related to the {311} 
defects and the formation and dissolution of the defects, from the interstitial clusters to the 
extended {311} defects, can be explained by common mechanisms. 

The proposition that the {311} defects are formed by condensation of interstitials and 
provide interstitial sources during boron TED has been supported by experiments using 
a variety of procedures to inject Si interstitials into bulk systems: the {311} defects are 
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observed to be formed by surface oxidation,^ and by GeV-electron irradiatiorJ13 and by ion 
implantations!! 

The {311} defects are often called rod-like, because they are typically elongated along 
the (Oil) direction as much as a micron.0'0 The width of the {311} defects ranges 1 - 100 
nm along the (233) direction, perpendicular to the elongation direction. The name {311} 
defects indicates the observed habit planed on which the rod-like defects lie, namely the 
{311} plane formed by the (Oil) and (233) directions. 

Here, we present total energy calculations based on the tight-binding Hamiltonian scheme 
(i) to study the structural properties and energetics of the {311} defects as function of their 
dimensions and interstitial concentrations and (ii) to find possible mechanisms of interstitial 
capture by and release from the {311} defects. We investigate the {311} defects system- 
atically from few-interstitial clusters to planar defects. Our results can be summarized as 
follows. 

(1) Interstitial chain - the basic building block of the {311} defects: we show that an 
interstitial chain along the (Oil) direction is stable with respect to isolated interstitials. 
The formation energy of the interstitial chains, E{ nt = 2 eV per interstitial, is smaller than 
the formation energies of isolated interstitials, 3-5 eV. Interstitial chains constitute basic 
building blocks of the {311} defects, i.e., extended defects on the {311} habit plane can 
be constructed by arranging the interstitial chains along the (233) direction (Sec. Ill and 
IV. A). 

(2) Stability of the {311} defects: interstitial chains along the (Oil) direction are stable 
against the isolated interstitials if they contain more than 2 interstitials. The formation 
energy per interstitial decreases linearly with the length of the interstitial chains along the 
(Oil) directions. More stable extended {311} defects than isolated interstitial chains are 
formed by condensation of interstitial chains along the (233) direction. The interstitial 
concentration of the most stable {311} defect, a planar defect, is 5 x 10 14 /cm 2 (Sec. Ill, 
IV. B and VII). 

(3) Growth mechanism: the stability dependence of the {311} defects on the interstitial 
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concentration indicates that finite-size interstitial clusters capture interstitials to grow into 
the interstitials chains along the (Oil) direction. The elongated interstitial chains are then 
further stabilized by capturing interstitials or interstitial chains side by side along the (233) 
direction. This growth mechanism based on the stability study is in a good agreement with 
the experimental claims that the elongation of the rod-like {311} defects along the (Oil) 
direction precedes the growth in the width along the (233) direction (Sec. IV.B). 

(4) Propagation of the interstitial chains: we propose a mechanism which can account 
for the motion of interstitial chains in the direction perpendicular to the chain direction 
with a relatively small energy barrier. Successive rotations of pair atoms on the {011} plane 
displace the interstitial chains. The growth of the {311} defects along the (233) direction can 
be explained by propagations of interstitial chains which are attracted to and captured by 
the already existing {311} defects. We show that V-shaped bend structure and atomic steps 
found in the transmission electron microscopy (TED) images can be formed by combinations 
of the interstitial chains and the planar rotations (Sec. V and VI). 



II. CALCULATIONAL DETAILS 

We perform total energy calculations based on a linear-scaling method in a tight-binding 
representation,^! using supercell methods at the T point. A tight-binding Hamiltonian de- 
veloped by Kwon et a/.@ is used to study the defect structures. This TB Hamiltonian gives 
a good description of the relative energies and equilibrium volumes of the diamond struc- 
ture and the metallic phases.0 More importantly, this TB Hamiltonian describes the elastic 
properties of the diamond structure with an error less than 5% compared to experiments 
and gives the formation energies of the point defects such as vacancies and interstitials which 
are in good agreements with those by the LDA calculations-^^ The validity and accuracy 
of this TB Hamiltonian to describe the diamond structure has been addressed by several 
calculations which show good agreements with the LDA calculations: reconstructions of the 
Si(100) surface^; the 90° partial dislocations in Sii and hydrogenated amorphous Si£3 
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The total energy calculations of defects requires large supercells to obtain converged 
formation energies, mainly due to the long-range structural relaxation. We use the orbital- 
based linear- scaling method implemented on Cray-T3D.il We choose a spherical localization 
of 6 A for less than 1% error in the total energy (about 30 meV per atom) compared to 
those obtained by exact calculations by diagonalizations. The chemical potential, /i, is 
adjusted to achieve the correct number of electrons with an error less than 10 -5 electron 
charge. Significant charge transfer is prevented by using a finite Hubbard-like term ({7 = 4 
eV); however, little differences are observed from a finite U and U = in the relaxed atomic 
structures and in the total energies. 

Figure 1 shows the smallest orthorombic unit cell used to study {311}defects. The [311] 
direction, normal to the habit plane of the defects, is chosen as the z axis and the [Oil] and 
[233] directions as the x and y axes, respectively. The unit lengths along the three principal 
axies are L xo = a/y/2, L yo = a\/Tl/\/2 and L zo = a\/TT, for the lattice constant a of the 
diamond structure Si. The lengths of a computational cell along the [Olljand [233] directions 
are varied by choosing integer multiples (n x ,n y ) of L xo and L yo and the length along the 
[311] direction at L z = 2L zo . Values for n x and n y are chosen so that the displacement of 
the atoms far from the defect core is less than 0.02 A with respect to the regular lattice sites 
of the perfect diamond structure. Periodic boundary conditions are applied along all three 
directions. For the structural optimizations, initial configurations of the model structures 
are given properly and constant temperature MD simulations are performed at 300 - 600 
K for about 1 psec to avoid trapping in the local minima of the atomic structures with 
small energy barriers. Then, atomic positions are fully relaxed by using the steepest descent 
method until the atomic force on each atom is less than 0.01 eV/A. The effective temperature 
of the relaxed ionic configuration is less than 0.1 K. 

We define formation energy of a defect structure which contains N int interstitials and 
Nbuik Si atoms in a computational cell as 

E f = E tot [N int + N bulk ] — E t ot[Nbuik], (1) 

^bulk 
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The formation energy per interstitial, E{ nt = /N int , is used to compare defects at a wide 
range of interstitial concentrations. A smaller E( nt corresponds to a more stable interstitial 
defect. Calculations using and using the TB HamiltonianEl found that the (110)- 

interstitialcy is the most stable point defect with a formation energy of 3.2 — 3.9 eV. Our 
total energy calculation using a 1000-Si supercell gives i?/ 110 \ = 3.9 eV. 

A binding energy of interstitials with respect to isolated interstitials is defined as 

- E b = E f - N int E( no) , (2) 

for Et 10 K, the formation energy of an isolated (llO)-interstitialcy. A positive binding energy 
of a defect containing interstitials indicates that the interstitial defect is stable with respect 
to isolated interstitials. 



III. INTERSTITIAL CHAINS ALONG THE (Oil) DIRECTION 

Here, we discuss an interstitial chain structure which is stable against isolated interstitials 
and which constitutes a building block of the {311}defects. The interstitial chain structure 
in Fig. 2 shows that interstitial chains along the (Oil) direction can be inserted into bulk 
Si without introducing any dangling bond by stacking pair interstitials with a periodicity of 
L xo . This arrangement of interstitials along the (Oil) direction is favorable, since only two 
dangling bonds would be required for any finite-length interstitial chainBH Since the number 
of dangling bonds per interstitial is inversely proportional to the length of the chain, the 
interstitial chain becomes more stable as it grows along the (Oil) direction. The interstitial 
chain structures have also been used for an atomic model of planar {311} defects by Takeda 
et al^E 

The most simple defect configuration containing interstitial chains, the structure in 
Fig. 2(a), is obtained by inserting an interstitial chain into bulk Si. By adding two intersti- 
tials per plane, two bonds (dotted line) are broken and two seven-members are introduced 
along the ±[311] directions. The six-member ring at the center turns into two adjacent 
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five-member rings. The bond-angle distortion ranges from -23° to 20° and the bond-length 
distortion from -0.13 to 0.03 A. This chain configuration is a stable interstitial complex 
with a relatively small formation energy per interstitial, E( nt = 2.2 eV, compared to that of 
isolated interstitials (3-5 eV). 

More stable structures can be obtained by eliminating five-member rings which share 
common bonds. A rotation of the atoms, connected by a solid bond indicated by an arrow, 
converts the right five- member ring in Fig. 2(a) into a six- member ring in Fig. 2(b) and lowers 
the formation energy to 1.7 eV. This defect structure is characterized by a six-member ring 
surrounded by five- and seven-member rings and has a mirror symmetry with respect to 
the center of the six-member ring.ll While the concentration of the additional atoms is two 
per unit length (L xo ), this configuration appears to have two interstitial chains due to its 
symmetry. 

A rotation of the other solid bond results in an interstitial chain surrounded by six- and 
five-member rings in Fig. 2(c). This configuration has an inversion symmetry with respect 
to the center of the interstitial pairs. The formation energy is 1.7 eV per interstitial and 
the corresponding bond-angle and bond-length distortions are —15° to 20° and -0.1 to 0.03 
A. We label the structure in Fig. 2(b) as an /-chain on {100} plane and that in Fig. 2(c) 
as an /-chain on {311} plane, based on the symmetry of the interstitial chains. These 
stable interstitial chains constitute basic building blocks of the {311}defects. Arranging 
interstitial chains side by side in the (233) direction results in the formation of even more 
stable extended defects lying on the {311}plane. 

The small formation energy per interstitial (< 2 eV) of the interstitial chain structures 
in Figures 2(b) and (c) suggests that microscopic defects containing only a few interstitials 
may exist as stable structures as well. Interstitial clusters containing 2-6 interstitials are 
constructed by inserting finite-length interstitial chains into bulk Si. Table I gives the 
formation energies per interstitial (E{ nt ) of finite-size interstitial clusters. Generally, E{ nt 
decreases as the size of interstitial defects increases, while the binding energy increases. 
This E int dependence on the number of interstitials, i.e., the length of the interstitial chain, 
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is consistent with experimental observation that the rod-like defects are extended along the 
(Oil) direction as long submicron. EMS The elongation of the {311}defects in the (Oil) 
direction is the consequence of the formation of the energetically favorable interstitial chain 
structures. This observation agrees with the energetic argument based on the minimum 
dangling bond ratio of the (Oil) chain structures. 

Table I shows that microscopic defects containing more than two interstitials are stable 
with respect to isolated interstitials. A binding energy per interstitial can be interpreted 
as an average energy required to release interstitials from the interstitial clusters and the 
interstitial chains. The energy cost to evaporate interstitials increases as the size of the 
interstitial cluster grows: a smaller interstitial cluster would be dissolved at a lower annealing 
temperature or in a shorter period of time. Recent experiments reported boron TED with 
a small activation energy (< 2 eV) even when no macroscopic defects are detected. 1. The 
lack of extended defects in their samples was attributed to the low-energy ion-implantation 
condition. It is possible that few-interstitial clusters or interstitial chains, which are stable 
with respect to isolated interstitials, are generated during the low-energy ion implantation 
and contribute to the boron TED by providing interstitials during the thermal annealing. 



IV. EXTENDED ROD-LIKE {311}DEFECTS 
A. Structure 

Here we show that extended {311}defects can be constructed by arranging the interstitial 
chains along the (233) direction. It has been observed that the rod-like defects elongated 
along the (Oil) direction grow thicker along the (233) direction, when there is a constant 
supply of interstitials.lll Possible configurations consisting of two interstitial chains are shown 
in Fig. 3. Extended {311}defects can be constructed by adding more interstitial chains along 
the [233] direction. As seen in Fig. 3, the {311}habit plane of the defects is made of two 
orthogonal directions, the (Oil) and (233) directions. In forming defects extended along 
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the [233] direction, the arrangement of interstitial chains should satisfy following conditions: 
when two interstitial chains are separated by L yo /2 as in Fig. 3(a), the pair include an I* and 
an I*\ when separated by L yo , the pair consist of two identical interstitial chains, two I*'s as 
in Fig. 3(b) or two /J-'s. Otherwise, dangling bonds are introduced and the corresponding 
distortions in bond angles become much larger than typical bond-angle distortions of ±20° 
of the stable {311}defects. Indeed, molecular dynamics simulations at 600 K about 1 psec 
remove such configurations. 

A new structure, an eight-member ring, is shown in Fig. 3(b). When two interstitial 
chains are separated by L yo in the [233] direction, there are three possible configurations: 
an eight-member ring (Os); a seven- member ring (Or); and a six- member ring (Oq). The 
three structures between two interstitial chains are related to each other by the same kind 
of planar rotations in Fig. 2. The eight-member ring in Fig. 3(b) can be converted into 
an O7 or an Oq unit by rotations of the solid bonds (heavy lines) on the {011} plane. An 
atomic model proposed by Takeda uses the eight-member ring as a basic unit in constructing 
the planar defects.ilSS However, we consider three possibilities, the Os, O7 and Oq units 
between two interstitial chains separated by L yo . For all the structures we have studied, the 
six-member rings (0§) between interstitial chains are found to be unstable with respect to 
the Os or O7 units. A structure which forms a boundary between the defects and the bulk 
Si is denoted as E 7 according to the notation, introduced by TakedailE Similarly to the 
eight-member rings, the seven-member rings (E?) can be converted into six-member rings 
(Eq) by planar rotations of the solid bonds in Fig. 3. 



B. Energetics 

The formation energies of rod-like {311}defects containing few interstitial chains are 
listed in Table II. For example, the E^IIOgHE? model structure has the elements: / - 
interstitial chains; Os - eight-member rings between the interstitial chains; and Ej - seven- 
member rings at the defect boundaries. Among possible combinations of the O and E units 
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for the given arrangement of interstitial chains, configurations with a lower formation energy 
are listed in Table II. The rod-like defects become more stable with increasing number of 
interstitial chains (I units) in the [233] direction as seen the E^s of Table II: the formation 
energy decreases (i) when interstitial chains are added side by side in the [233] direction 
with a distance L yo /2, Table II (a) — > (e); and (ii) when interstitial chains are added in the 
[233] direction in the presence of an O unit, Table II (f) — > (i). 

As seen in Table II, the eight-member rings (0 8 ) are stable when they are separated 
by a distance equal to or larger than L yo . The formation energy of the E 7 IO$IE 7 model 
structure is lowered by 0.1 eV by transformation of the 0$ unit to a seven-member ring (0 7 ). 
On the other hand, the eight-member ring of E 7 II0 8 IIE 7 model is stable with respect to 
the E 7 II0 7 IIE 7 structure. Similarly, we observe that the seven-member ring at a defect 
boundary is stable, if there are more than 2 / units inserted between E 7 and the closest O 
or E unit. For example, the total energy increases when the seven-member rings (-EV's) are 
converted into six-member rings for the E 7 IIE 7 and E 7 II0 8 IIE 7 models. 

We define an energy release AEj, when an interstitial chain is added to a rod-like defect 
X containing N int interstitials per L xo , as 

AEj = (N mt x EL [X] + 2 x E{ nt [E 7 IE 7 ] ) - (N int + 2) x E s mt [X + I]. (3) 

The formation energy, E- nt [E 7 IE 7 } = 1.8 per interstitial, is used as the formation energy 
of an isolated interstitial chain, which contains two interstitials per unit length L xo . The 
energy release AEj for the E 7 IIO%IE 7 structure to capture an interstitial chain and become 
the E 7 II0 8 IIE 7 structure is 1.6 eV per unit length (L xo ). On the other hand, A_E/ = 1.0 
eV for the E 7 IIIE 7 structure to become the E 7 IIIIE 7 structure. The larger AEj for the 
E 7 IIO$IE 7 structure indicates that the E 7 IIOgIE 7 structure provides a more efficient sink 
of interstitials than the E 7 IIIE 7 structure, which contains the same number of interstitials 
per unit cell. The A£/'s for the E 7 IIE 7 structure suggest that the interstitials can be 
evaporated with a smaller energy cost from the E 7 IIOgIE 7 structure than from the E 7 IIIE 7 
structure. 
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Growth mechanism of the {311}defects are suggested by both the E( nt at different inter- 
stitial concentrations and the energy release AEj , when an interstitial chain is added to an 
existing {311} defect. Previously we showed that an interstitial chain can exist as a stable 
structure against isolated interstitials with an E{ nt < 2.2 eV, which is smaller than that of 
isolated interstitials. The interstitial chain can provide a sink of interstitials to capture in- 
terstitial chains along the [233] direction. Reactions which release a larger energy to capture 
an interstitial chain are more likely to happen during the growth process. 

We express the capturing of an interstitial chain by a rod-like {311}defect symbolically 



as 



■■■IE + I — ► • • ■ IIE 7 for d = L yo /2, 

■■■IE + I — ► • ■ ■ IOIE for d = L yo . (4) 

Here, d is the distance of the additional interstitial chain from the boundary interstitial chain. 
Table II shows that the energy release by addition of an interstitial chain is similar for either 
d = Ly /2 or d = L yo , which suggests that the interstitial-chain capturing mechanisms in 
Eq. (|D are equally likely to happen. In other words, an interstitial chain can be added next 
to the boundary interstitial chain with a distance either L yo /2 or L yo . Consequently, the O 
units would be randomly introduced and the {311}defects would have no particular periodic 
arrangement of the interstitial chains. In fact, experiments observe no particular periodicity 
of the {311}defects in the (233) direction,!! while the periodicity in the (Oil) direction is 
identified to be L xo , which results from the interstitial chain structure. 

We denote the "interstitial density" to be ratio of the number of interstitial chains to 
the total number of units in the defects: 

"interstitial density" = _ - ^ - ^ (5) 

We obtain the optimal ratio of the interstitial chains of 67% for the {311}defects based on 
following results: (i) the formation energy for the extended defects (3) - (5) in Table II 
is same with differences less than 0.1 eV. This formation energy, E int = 1.3 eV, is smaller 
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than E int = 1.4 eV of a planar defect - an infinite sequence of the interstitial chains in the 
[233] direction (■ • -IIII ■ ■ ■ = /I/), (ii) The planar defect consisting of sequencies of 110$ 
structure, • • ■ II0 8 II0 8 • • • = /II0 8 / model, is found to be very stable with the smallest 
Ef nt = 1.2 eV among the model structures we have investigated. This interstitial density 
(67%) is in a good agreement with the observed ratio of / units, 62% from HRTEM images.il 



V. PROPAGATION OF INTERSTITIAL CHAINS 

So far we have discussed the energetics of classes of defects, from the interstitial clusters 
to the extended {311}defects, solely in terms of the formation energies. It has been shown 
that the interstitial clusters can bind interstitials and become elongated interstitial chains 
along the (Oil) direction. The interstitial chains in turn grow to the extended {311}defects 
by capturing interstitials along the (233) direction. We propose that successive planar 
rotations, which introduce eight-member rings between interstitial chains separated by L yo , 
can also displace interstitial chains. By successive rotations of atoms connected by solid 
bonds in Fig. 4, an interstitial chain characterized by the {311}habit plane is displaced by 
A = 2a/v / 2 along the arrow (the [Oil] direction), perpendicular to the [Oil] chain direction. 
The planar rotations connect two interstitial chain structures with the same E( nt = 1.7 
eV - the /-chain on the {311}plane, Fig. 4(a)(c)(d), and the /-chain on the {100} plane, 
Fig. 4(b)(e). 

This kind of coordinated atomic motions was first introduced by Pandey as a diffusion 
mechanism of Si atoms in a thermal equilibrium without introducing point defects.il He 
suggested that successive rotations of nearest neighboring atoms in the {011} plane and the 
{011} plane result in exchanges of atoms and eventually displacements of Si atoms. The 
interstitial chain configurations shown in Fig. 4 are related to each other by (f> ~ 70° rotation 
of the pair atoms (solid circles) in the {011} plane. 

For E* , a formation energy of the interstitial chain, Eb, an energy barrier for the prop- 
agation, and S, an entropy involved with the propagation, the probability for an intersti- 
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tial chain to propagate is proportional to exp { — (E- f + E b — TS)/k B T}. The propagation 
mechanism of the interstitial chains via planar rotations has such desirable features: (1) no 
dangling bond is introduced during the transformation due to atomic relaxations and thus 
a small energy barrier is expected; and (2) the entropy involved with the rotation may be 
high because of the large extent of atomic relaxations and many possible paths. 

An energy barrier for the transformation from the structure of Fig. 4(a) to (b) is estimated 
in two ways. First, an energy barrier for the rotation is estimated by calculating total 
energies of saddle points and comparing the total energies with those of the interstitial chain 
structures. The calculations are performed at a fixed computation cell of 4L xo x 3L yo x 2L zo . 
This computational cell contains four solid bonds into the {011} plane which are not shown 
in Fig. 4(a). The saddle point configurations are initially given by simultaneously rotating 
the atoms connected by the solid bonds by <fi = 15°, 30°, and 45° on the {011} plane. The 
bond length of the two atoms is chosen at the average value of the bond lengths of two 
structures in Fig. 4(a) and (b). Then, structural relaxations of the initial configurations are 
performed by allowing all atoms to move except for the two rotated atoms. The rotated 
atoms are fixed at the intermediate positions. The calculated energy barrier is 1.7 eV per 
bond. When the similar rotation of two atoms is performed in the bulk Si without any 
interstitial, the energy barrier is about 2.7 eV per bond. 

Alternatively, the energy barrier is estimated by rotating one bond at a time among four 
solid bonds in a computational cell, which are underneath of the solid bond in Fig. 4(a) 
in the [Oil] direction. For rotations of one to three bonds, the total energy increases by 
1.2 ± 0.1 eV in comparison with those of the interstitial chains in Fig. 4(a) and (b). The 
corresponding increase in the formation energies is less than 0.2 eV per interstitial. No 
dangling bond is introduced by these planar rotations and the intermediate structures have 
bond-angle distortions of ±22°, slightly larger than —15° to 20° bond- angle distortions of 
the interstitial chain configurations in Fig. 4(a) and (b). 

The energy barriers estimated by two different paths suggest that many possible combi- 
nations of the planar rotations, which has an energy barrier less than 2 eV per bond, can 
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lead to the displacements of interstitial chains. The growth of the {311}defects along the 
(233) direction during thermal annealing0 can be explained by the coalescence of interstitial 
chains which propagate via the planar rotations. The same mechanism may be responsi- 
ble for the unfaulting of the {311}defects which relates the {311}defects with the perfect 
dislocation loops at high temperatures (> 1000 K). ? 

VI. STEP STRUCTURES AND V-SHAPED BENDS 

We now show that step structures and V-shaped bends found in the extended 
{311}defects can be constructed by combination of rotations of nearest neighboring atoms 
on the {011} plane together with the interstitial chain structures. Rotations of three solid 
bonds in Fig. 5(a) result in a displacement of the right interstitial chain. The formation en- 
ergy of the structures before and after the transformation is the same at E{ nt = 1.5 eV. The 
structure in Fig. 5(b) is similar to the core model of the atomic steps observed in HRTEM 
images.ii Addition of interstitial chains to this defect structure along ± y direction (the 
[233] direction) leads to the formation of a defect structure containing an atomic step. The 
two {311}habit planes shown as dashed lines are separated by a step height of A = aj\f\\ 
along the {311}direction, which agrees with that predicted in Ref.[25]. 

It has been observed that the {311}defects change habit planes, for example, between 
the {311} and {311} planes, i! The atomic configuration around the bend where the habit 
planes change can be identified as the /-chain in the {100} plane in Fig. 2(b). We show in 
Fig. 6 that the interstitial chains can be added to an /-chain {100} plane to form a V^-shaped 
defect. The initial configuration, when two interstitial chains are added, is chosen so as to 
avoid the generation of adjacent five-member rings by bond rotations on the {011} plane. 
The lowest energy configuration of the bend structure is found to have seven-member rings 
at the edges {E 7 ): by transformation of the six-member rings to E 7 J s, the total energy is 
lowered by 3.8 eV per cellP corresponding to AE{ nt = —0.2 eV. The bond-angle distortion 
ranges between -17 ° and 26 ° and the bond-length distortion between -0.1 A and 0.02 A. 
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The formation energy, 1.3 eV per interstitial, is comparable to that of the EjIIIEj model 
structure at the same interstitial concentration of 6 per unit length, L xo = a/\^2. The 
defect can grow on both habit planes, the {311}and {311} planes, by capturing interstitials 
or interstitial chains. The stability of the E?s suggests that the additional interstitial chains 
are to be separated by L yo . In other words, for the center interstitial chain /, the V-shaped 
vend structure is likely to have an interstitial arrangement as • • • IO%IIIO% ■ ■ ■. 

VII. PLANAR {311}DEFECTS 

Extended {311}defects as wide as 100 nm along the (233) direction are studied by ap- 
proximating them as planar defects which are periodic both in the (Oil) and the (233) 
directions.il Formation energies of planar defects, constructed by combinations of intersti- 
tial chains (/) and eight-member rings {Os) along the [233] direction, are listed in Table III. 
As introduced by Takeda et a/.,0 the periodicity of a planar defect along the [233] direction 
is specified by an arrangement of J's and O's within a unit cell, which is denoted by / • ■ ■ /. 
For example, the jllO^j model contains two / units separated by a distance of L yo /2 in 
the [233] direction and one 0§ unit within a unit cell. The periodicity of the jllO^j model 
along the [233] direction is 3(L yo /2), where L yo /2 is the average width of each unit and the 
factor 3 is the total number of units within a unit cell. 

The formation energy per interstitial of the planar {311}defects is significantly smaller 
than that of isolated interstitials and comparable to that of the extended {311}defects 
with finite widths along the [233] direction (Table II). The stability of the planar defects 
against isolate interstitials indicates that the {311}defects can grow indefinitely as long as 
the interstitials are supplied, e.g., during the ion implantation process. Among the model 
structures we have investigated, the jllO^j model is the lowest in the formation energy per 
interstitial (E( nt = 1.2 eV). The density of the / units of the jllO^j planar defect is 67%, 
corresponding the interstitial concentration of 5 x 10 14 /cm 2 = 4/(L xo x 3(L yo /2)). 

The model structure with adjacent 8 units along the [233] direction is energetically 
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least favorable at the same interstitial concentration in comparing the /IlOg/, /IHOglOs/ 
and / IIIIO$Os/ models. In fact, segments including two adjacent eight-member rings, 
■ ■ -OsOs ■ ■ •, have not been identified as local structures of the extended {311}defects.@il 
When a planar rotation is applied on the common bond shared by the neighboring eight- 
member rings of the jllllO^O^j model, E int is lowered by 0.2 eV. The stability of the 
eight-member rings of each structure is studied by applying planar rotations of the atoms 
which form boundaries between the eight-member rings and six-member rings along the 
[233] direction. The last column of Table III shows that the O^s are stable when their 
separation is equal to or larger than L yo , consistently with the result of the finite-width 
{311}defects. 

The difference of the E{ nt 's of the /IlOg/ and /IHOglOg/ models implies that the 
stability of the planar defects is determined not only by the density of interstitial chains but 
also the arrangement of interstitial chains. Planar defects containing isolated interstitial 
chains between O units (• ■ ■ OIO ■ ■ ■) have higher formation energies than those without 
them at the same interstitial concentration. The energy gain by transformation of one of 
the eight-member rings to a seven-member ring is insufficient to change the relative stabilities 
of the planar defects. Based on the E- nt of the planar defects listed in Table III, we expect 
that the / IIIO%IIO%/ and jlllO^j planar defects to have comparable E int to the most 
stable planar defect, the /IIO$/ model. 

The formation energy of the planar defects we have discussed is obtained without al- 
lowing relaxation of the supercells. A supercell relaxation would have generated the dis- 
placement vector (Burger vector) around the core of the {311}defects identified by TEM 
observations.0'@ However, we anticipate that the relative stability is correctly described by 
the fixed cell calculations. The critical dimension to obtain the converged formation energy 
for the planar defects is L z , the length of the supercell in the [311] direction. We obtain the 
formation energies of the planar defects at L z = 3L zo and the formation energies are also 
listed in Table III. The formation energy is lowered by 0.08 eV on average for all the planar 
defects. The relative stability of the planar defects is the same for L z = 2L zo and 3L zo . 
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In Fig. 7, we show the formation energy of the defects which are infinite along the [233] 
direction and finite (L D = n D L xo ) along the [Oil] direction. The elongation direction of the 
{311}defects is always along the (Oil) direction. However, we study the hypothetical defect 
structures with lengths (L D ) shorter than the width in the [233] direction to compare with 
the planar defects. Since the infinite limits of these defects (Lb oc oo) are planar defects, 
we can compare the formation energy of the planar defects in Table III with E( nt [oo], which 
can be extracted from the finite-length-defect calculations. 

The defect structures in Fig. 7 are constructed by inserting a finite-length planar defect 
between bulk layers. The defects become either the /// or jIO%j model structure at the 
infinite limits. The formation energy increase linearly to the length along the [Oil] direction: 

Ef = E f edge + N mt EUoo\. (6) 

Then the E- nt of the defects with respect to the length is written as 

E{ nt = eUoo] + E{ dge /N int oc 1/L D (7) 

The slopes, E{ nt [oo], of two kinds of defects is the formation energy per interstitial at the 
infinite limits, i.e., N int oc L D — > oo. The slopes obtained from linear fittings are in good 
agreements with the E{ nt 's in Table III, 1.4 eV and 1.7 eV for the /// model and /IO$/ 
model, respectively. This result confirms that the /// model is more stable than the /I0 8 / 
model. 

Finally, the ability of the {311}defects to capture or release interstitials is studied by 
calculating an energy required to add an interstitial to the bulk Si in the presence of the 
planar defects (Fig. 8). An interstitial is added at the hexagonal site at a distance of h 
from the {311}habit plane. The eight-member rings of the jlO^j planar defect can provide 
efficient sinks of interstitials: the extra energy to add an interstitial to the /I0 8 / planar 
defect decreases with decreasing h, 3.2 eV at 16.5 A and 1.5 eV at 10.1 A. When the 
interstitial is added right above the eight-member ring (3.8 A from the habit plane), about 
0.1 eV energy is released. On the other hand, an additional interstitial is not bounded by 
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the /IlOg/ planar defect. The energy required to add an interstitial is 4.7 - 5.1 eV. This 
energy is comparable to the formation energy of an isolated interstitial at the hexagonal 
sites and larger than that of the (110) interstitialcy. These results support the calculations 
showing that the saturation of interstitials is achieved at the / density about 67% in the 
/IIO a / structure (Table III). 

VIII. CONCLUSIONS 

We have showed that interstitial chain structures elongated along the (Oil) direction 
constitute the basic building blocks of the extended {311}defects. This interstitial chain 
configuration is favored due to the minimum ratio of dangling bonds per interstitial and 
small distortions in bond angles and bond lengths. Even finite-size interstitial chains, as 
small as 2-interstitial clusters, are found to be stable against isolated interstitials. A growth 
mechanism of the {311}defects can be speculated from the dependence of the E^s on the 
dimensions of the defects and the interstitial concentrations. The growth of the {311}defects 
would occur first by (1) elongation of interstitial chains along the (Oil) direction and (2) the 
widening of the defects follows by capturing interstitial chains alongside the (233) direction. 
We present an efficient mechanism, rotations of atoms on the {011} plane, which can lead 
to propagation of interstitial chains. The planar rotation together with the interstitial 
chain structures can be applied to explain growth process of the {311}defects and related 
structures such as the V-shape bend structure and atomic steps. We observed that the most 
stable arrangement of interstitial chains of the extended {311}defects is ■ ■ ■ II0 8 II0 8 - ■ ■ 
and the optimal density of interstitial chains is 67% which is in a good agreement with the 
experimental value of 62%. 
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FIGURES 



Fig. 1 An orthorombic unit cell used to describe the {311} defects. The three principal 
axes are along the [Oil], [233] and [311] directions. For the Si lattice constant (a = 5.43 
A), the dimensions of the unit cell are L xo = a/y/2,L yo = aVTl/y/2 and L zo = ay/Tl. 
Supercell calculations at the T point are performed by taking a computational cell whose 
size is L x = ti x x L xo , Ly — Tiy x L yo , and L z = 2 x L zo . Periodic boundary conditions are 
applied in all three directions. 

Fig. 2 An interstitial chain obtained by stacking pair interstitials with a periodicity of 
L xo along the (Oil) direction. The solid atoms are L xo /2 into the plane with respect to 
the open-circled atoms, (a) An interstitial chain inserted into bulk Si is surrounded by two 
adjacent five-member rings. The dotted lines denote the bonds which are broken by addition 
of interstitial chains. This structure has a formation energy of 2.2 eV per interstitial, (b) By 
a rotation of atoms connected by a solid bond (arrow), a more stable structure with E( nt = 
1.7 eV is obtained. The structure has a mirror symmetry with respect to the center of the 
six-member ring. The dashed line indicates the {100} plane normal to the center line, (c) 
A rotation of the other bond results in an interstitial chain configuration of E{ nt = 1.7 eV. 
The symmetry of the structure is an inversion symmetry with respect to the center of pair 
interstitials on the {311} plane (dashed line). 
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Fig. 3 Projections on the {011} plane of the rod-like {311}defects containing two inter- 
stitial chains. The solid atoms (•) are deeper into the plane than the open-circled atoms 
(o) by L xo /2. Interstitial chains are specified by the upper atoms with respect to the {311} 
plane (dashed lines): (1) l\ when the upper atoms are out-of the {011} plane; and (2) l\ 
when the upper atoms are into the {011} plane. The rod-like {311}defects are believed to 
grow by adding interstitial chains along the [233] direction. Total energy calculations give 
the formation energy per interstitial at 1.5 eV for (a) and at 1.6 eV for (b). Seven-member 
rings (-EV's) form boundaries between the defects and the bulk Si, and they can be converted 
into six-member rings (E 6 ) by rotations of the solid bonds (heavy lines). In a similar way, 
an eight-member ring 8 can be converted into either a seven-member ring (0 7 ) or a six- 
member ring (0 6 ). The stability of the E 7 and 8 units depends on the arrangement of 
interstitial chains. A rotation of any solid bond of (a) increases the total energy. However, 
more stable structures can be obtained by rotations of the solid bonds in (b) as far as the 
transformation doesn't introduce adjacent five-member rings. 

Fig. 4 By succesive rotations of solid atoms connected by solid bonds, an interstitial chain 
on the {311} habit plane is displaced by A = 2aj\[2 along the perpendicular direction of 
the interstitial chain (arrow). The <g> denotes a reference atom. Intermediate structures are 
characterized by the habit planes, the {100} plane for (b) and (d) and the {311} plane for (c), 
indicated by dashed lines. The formation energy per interstitial of the illustrated structures 
is 1.7 eV for a fixed computational cell. The dimension of the cell is AL xo x 3L yo x 2L zo . 
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Fig. 5 By rotations of three solid bonds, a rod-like {311}defect containing two /-chains 
on the {311}plane (a) is transformed to a step configuration with a height of A = a/y/Ti 
along the [311] direction (b). The solid atoms are deeper into the plane by L xo /2 than the 
open-circled atoms. The difference of the total energies per unit cell of the two structures, 
which contains 528 bulk Si atoms, is less than 0.1 eV, and the corresponding AE( nt is less 
than 0.05 eV. Intermediate states after rotations of one or two solid bonds have the same 

jC 'int- 

Fig. 6 Bend configuration of the {311}defects (bottom) obtained by addition of two 
interstitial chains to an /-chain on the {100} plane (top). Insertion of interstitial chains 
results in breaking of bonds (dotted lines). The bottom structure contains 6 interstitials 
per unit length of L xo : 2 interstitials associated with the /-chain on the {100} plane and 
4 additional atoms indicated by arrows. A row of the center six-member rings along the 
[Oil] direction forms a boundary of the defects which have the {311} and {311} habit planes 
(dashed lines). This V-shaped defect has a mirror symmetry with respect to the center of 
the six-member rings. The formation energies per interstitial are 1.7 eV (top) and 1.3 eV 
(bottom). Further growth of the defect can occur by addition of interstitial chains on the 
both habit planes. 
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Fig. 7 Formation energy (E?) per unit width (L yo ) with respect to the number of 
interstitials per unit width (left); a schematic diagram to show how the computational cell 
is constructed (rig ht). A slab with a finite length of noL xo , for an integer no = 1- • - 4 
, contains defect structures based on the /I0 8 / and /// models. The slab is sandwiched 
between Si bulk layers along the [Oil] direction with a fixed total length of QL xo . A common 
width along the [233] direction is chosen at L yo . The number of interstitials per unit width 
(L yo ) is proportional to the length of the defects: N int = An D for the defects based on the /// 
model and N int = 2ud for those based on the /IO$/ model. The formation energy is fitted 
by Eq. (6). For the defects based on the /// model, the values of E{ dge and E( nt [oo] are 8.4 
and 1.4 eV respectively; E{ dge = 6.1 eV and E{ nt [oo] = 1.7 eV for the defects derived from 
the /IO s / model. The formation energy per interstitial decreases as the length increases as 
in Eq. (7). At N int — > oo, the E{ nt 's converge to those of the planar defects, 1.4 eV for the 
/// model and 1.7 eV for the /IO a / model. 

Fig. 8 Initial configurations to obtain the extra energy required to add one interstitial 
to the bulk system in the presence of a planar defect. The additional interstitial is located 
at a hexagonal site with a distance h measured from the {311} plane on which the center 
of interstitial chains lie (the dashed line). Two cases with hi = 16.5 A and h 2 = 3.8 A are 
illustrated at the same time. The computational cell is 6L xo x 3L yo x 2L zo . The distance of 
interstitials in the neighboring cells, i.e. the distances between images related by periodic 
boundary conditions, is larger than the distance of the additional interstitial to the habit 
plane of the planar defect. The interstitials associated with the interstitial chains of the 
/IOs/ model structure are indicated by arrows. 
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TABLES 





/-chain 


on {311} plane a 


/-chain 


on {100} pla 


ne b 


N int \ eV 


Ef/N int 


E b 


E b /N mt 




£ b 


E b /N int 


2 


4.7 


-1.6 


-0.8 


4.9 


-2.0 


-1.0 


3 


3.7 


0.6 


0.2 


3.4 


1.5 


0.5 


4 


3.4 


2.0 


0.5 


3.3 


2.4 


0.6 


5 


3.0 


4.5 


0.9 


2.9 


5.0 


1.0 


6 


2.8 


6.6 


1.1 


2.8 


6.6 


1.1 


oo 


1.7 




2.2 


1.7 




2.2 



Orientation of the computational cell: a [011] x [322] x [311] 



6 [011] x [011] x [100] 

TABLE I. Formation energy per interstitial, Ef nt in Eq. (1), and the binding energy, E b in 
Eq. (2), of few-interstitial clusters. The interstitial clusters are constructed by inserting interstitial 
chains between bulk layers whose length is at least 6 x L xo along the [011] direction. The length of 
an 7Vj n t-interstitial cluster along the (011) direction is Ni n tL xo /2. The orientation of the computa- 
tional cells is chosen according to the symmetry of the interstitial chains. The bottom row shows 
E( nt and E b /Ni nt for the line defects- infinite interstitial chains along the [110] direction (Fig. 2). 
The size dependence of E int 's shows that the interstitial clusters can trap interstitials to form more 
stable elongated structures. 
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Model structure 


N int (% of / units) 
per unit cell* 


E* 
eV 


AEj 


:X + / -► y 


(a) E 7 IES* 


2 (67) 


1.8t 


1.2: 


(a)+J -► (b) 








0.9: 


( a ) + / _^ (f) 


(b) E 7 IIE 7 


4 (50) 


1.5 


1.6: 


(b)+J -► (c) 








1.0: 


(b)+J -► (g) 


(c) E 7 IIIE 7 


6 (60) 


1.3 


1.2: 


( C )+J _> (d) 








0.6: 


( C )+J (h) 


(d) E 7 IIIIE 7 


8 (67) 


1.3 


1.0: 


(d)+J (e) 


(e) E 7 IIIIIE 7 


10 (71) 


1.3 






(f) E 7 I0 7 IE 7 or E 7 I0 8 IE 6 


4 (40) 


1.6 


1.3: 


m+7 ( R ) 


(g) E 7 II0 7 IE 7 or E 7 II0 8 IE 6 


6 (50) 


1.4 


1.1: 


(g)+/ - (h) 








1.3: 


(g)+/ - (i) 


(h) E 7 III0 7 IE 7 or E 7 III0 8 IE 6 


8(57) 


1.3 






(i) E 7 II0 8 IIE 7 


8 (57) 


1.3 







* Unit cell: L xo x QL yo x 2L zo 

** Identical to the /-chain on the {311}plane in Fig. 2(c). 
t = 1.7 eV for a unit cell L TO x 3L,,„ x 2L 20 



TABLE II. Formation energy of extended {311}defects consisting of interstitial chains. The 
computational cell is 3L xo x QL yo x 2L zo containing three unit cells in the [011] direction for all 
listed model structures. The ratio of the I unit to the extent of defects is given in percentage, 
e.g., 3/6 = 50 % for the E 7 II0 8 IE 7 structure. The energy release, AEj, is defined in Eq. 3 upon 
adding an interstitial chain (I) to the roddike defect X. Y denotes a defect structure constructed 
by adding an / with a distance of either L yo /2 or L yo from the boundary interstitial chain of the 
defect X. 
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Model structure 


mt 


EL(eV) 

int v / 


3L%o 


Stability of Os units 

J 


/I/ 


12 


1.35 


1.27 




110*1 


6 


1.68 


1.59 


Unstable, /I0 7 I0 7 I0 8 / 


/IIOsI 


8 


1.23 


1.16 


Stable 


/IIIOsIOs/ 


8 


1.35 


1.26 


Unstable, /III0 7 I0 8 / 


/IIIIIOs/ 


10 


1.24 


1.18 


Stable 


/IIIIOsOs/ 


8 


1.48 




Unstable, /IIII0 7 7 / 



* Per unit area: L xo x 3L. 



TABLE III. Formation energies per interstitial {E int ) of planar {311}defects consisting of in- 
terstitial chains and eight-member rings. The periodicity of the structures along the [233] direction 
is denoted by /• • •/. The number of interstitials and the formation energy are given per unit area, 
L xo x 3L yo , in order to compare structures with different periodicities along the [233] direction. The 
same computational cell is used for the total energy calculations of all the listed structures with 
two choices for the length of the cell in the [311] direction, 2L zo and 3L zo . The last column shows 
the stability of the eight-member rings and the most stable planar defects which can be obtained 
by rotations of atoms on the {Oil}. 
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Fig 1: Kim et al, Extended Si defects 




Fig 2: Kim et al, Extended Si defects 




Fig. 3: Kim et al, Extended Si Defects 




Fig. 4: Kim et al, Extended Si defects 



Fig 5: Kim et al, Extended Si defects 




Fig. 6: Kim et al, Extended Si defects 




Fig 7: Kim et al, Extended Si defects 




Fig. 8: Kim et al, Extended Si defects 



